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Abstract. - We study pressurised self-avoiding ring polymers in two dimensions using Monte 
Carlo simulations, scaling arguments and Flory-type theories, through models which generalise 
the model of Leibler, Singh and Fisher [Phys. Rev. Lett. 59, 1989 (1987)]. We demonstrate the 
existence of a thermodynamic phase transition at a non-zero scaled pressure p, where p = Np/4ir, 
with the number of monomers N — > oo and the pressure p — > 0, keeping p constant, in a class 
of such models. This transition is driven by bond energetics and can be either continuous or 
discontinuous. It can be interpreted as a shape transition in which the ring polymer takes the 
shape, above the critical pressure, of a regular N-gon whose sides scale smoothly with pressure, 
while staying unfaceted below this critical pressure. Away from these limits, we argue that the 
transition is replaced by a sharp crossover. The area, however, scales with N 2 for all positive p in 
all such models, consistent with earlier scaling theories. 



Fluid vesicles in three dimensions exhibit diverse shapes. 
This diversity originates in the nontrivial ground state 
configurations exhibited by a simple curvature Hamilto- 
nian for a closed surface as the enclosed volume and sur- 
face area are varied pQ. Not much is known about the 
behaviour of such a model system at non-zero tempera- 
ture. However, equivalent models in one lower dimension 
should provide useful insights into behaviour in the exper- 
imentally relevant three-dimensional case [UH]. The sim- 
plest such model is that of a self-avoiding polymer ring 
whose enclosed area couples to a pressure difference term 

Leibler, Singh and Fisher (LSF) [2] studied such a pres- 
surised polymer model, defining the ring in terms of teth- 
ered discs. The energetics of the model derived from a 
pressure term conjugated to the enclosed area. Monte 
Carlo simulations and a scaling analysis about zero pres- 
sure were used to demonstrate that the radius of gy- 
ration Rq and the averaged area A obeyed the scaling 
forms R 2 G ~ N 2u X{pN 2u ) and A ~ N 2v Y{pN 2l/ ), where 
v = 3/4 exactly in d = 2 and X(x) and Y(x) are scaling 
functions of a single scaling variable x. The dimension- 
less measure of the pressure difference is p — pa 2 /ksT, 
where p is the pressure difference, a is the diameter of the 
disc and T is the temperature, with the appropriate scal- 



ing variable defined as x — pN 3 / 2 . In a scaling regime 
around x ~ 0, data collapse for a range of system sizes 
was observed by LSF, with a x — transition separat- 
ing branched polymer behaviour, for p < 0, from a p > 
regime governed by the statistics of self-avoiding walks. 

The LSF scaling description applies for finite pN z / 2 . At 
large (positive) pressures, it is reasonable to expect that 
A ~ A^ 2 , with the expanded shape approaching a circle. 
However, extending LSF-type scaling arguments to the 
large pressure case yields A ~ A^ 3 for large pressures, as- 
suming only that pN <C 1 [B] . Any expanded regime with 
A ~ A^ 2 is thus inaccessible through an expansion about 
the LSF fixed point at p — 0. Related scaling arguments, 
based on a blob picture, have recently been presented for 
the case in which the scaled pressure p — Np/Air is finite 
and non-zero for large N [6j|7|. These suggest that self- 
avoiding pressurised chains should swell smoothly as p is 
increased, ruling out a phase transition at nonzero p. 

However, there are models for pressurised polymer rings 
which do exhibit a continuous pressurisation transition 
at non-zero p, (although these models allow for self- 
intersections) [9l[T0]. Recent work on ring polymers with 
a pressure term conjugated to the algebraic (signed) area, 
obtains a continuous transition, between a collapsed phase 
with (A) ~ A" and an expanded phase with (A) ~ A^ 2 
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Fig. 1: Configurations of a self-avoiding pressurised polymer, 
with finite beads connected by quadratic springs of maximum 
extension Ro (Model C). These figures illustrate sample con- 
figurations on either side of the transition; both for (a) p < 1 
as well as for (b) p > 1, where p is the critical pressure. 



[TT1[T2"] . The critical pressure, in a Flory-like theory, arises 
from the competition between an entropic R 2 /N term, 
which favours more compact configurations, and a pres- 
sure term —pR 2 /N, which favours expanded configura- 
tions. Such terms should also be present in theories for 
the self-avoiding case, supplemented by additional terms 
accounting for self- avoidance. Whether it is the true or the 
algebraic area which enters the Hamiltonian should also be 
irrelevant in the large pressure limit Thus, it is nat- 
ural to ask whether the inclusion of self-avoidance suffices 
to destroy the non-zero p transition altogether or whether 
some trace of the singular behaviour at p ~ 1 in the self- 
intersecting case survives. If it does, this would indicate 
a hitherto unanticipated feature of such pressurised poly- 
mer ring models, with possible applications to the three- 
dimensional case. 

This Letter demonstrates the existence of an unusual 
and somewhat subtle phase transition as p is increased 
across a critical value, in some models for pressurised, 
self-avoiding polymer rings. This transition can be either 
discontinuous or continuous, and separates a weakly in- 
flated phase from a strongly inflated phase, in a sense we 
make precise below. The average area (A) scales as N 2 on 
both sides of the transition, in accord with earlier scaling 
arguments. Crucially, the nature of fluctuations changes 
across the transition point, with the ring expanding with- 
out feeling the effect of the maximum extension Rq be- 
low the transition, while being constrained by Rq at and 
above the critical pressure. Sample configurations below 
and above the transition are shown in Fig. [TJ Interest- 
ingly, for models of the LSF type, such transitions appear 
to be absent in the thermodynamic limit, although their 
effects may be manifest in sharp crossovers. 

We first summarise prior scaling arguments for this 
problem [6, 7 . These use the blob picture, describing the 
polymer ring as a one-dimensional object fluctuating in 
two dimensions. In response to a relevant perturbation, 
the ring can be naturally decomposed into N/g blobs of 
g monomers each. A blob stores an average tensile en- 



ergy equal to fc^T = 1, i.e. 7^ ~ 1, where 7 is the in- 
duced surface tension due to the perturbation and £ is the 
blob size. At length scales smaller than the blob size £, 
the ring is unaffected by the perturbation, thus obeying 
£ ~ g v , where v is the swelling exponent characterising 
the unperturbed manifold statistics. This implies that the 
surface tension scales as 7 ~ g~ v ■ At length scales larger 
than £, the perturbation stretches the manifold; the to- 
tal projected length of the chain is the product of the 
number of blobs and the projected length of each blob: 
R ~ £ ~ A7 1 / l/ ~ 1 . For pressurised rings, the sur- 
face tension and the pressure difference p are related via 
Laplace's Law, *y/R ~ p. The total projected length then 

scales as R ~ jV^-ip 2 "- 1 , while the average area scales 
as 

(A) ~ R 2 ~ N 2 {pN)^, (1) 

with the surface tension given by 7 ~ pR ^ (pN) 2 "- 1 . 

The above scaling argument, while justified for pN <C 1, 
must be carefully examined for pN ~ I , i.e. when the blob 
size and the microscopic cutoff are comparable. For the 
blob picture to be valid, the number of blobs should be 
larger than unity as well as less than N\ these conditions 
impose N~ 2v < p < N^ 1 , or equivalently pN < I. 

The fact that such an analysis predicts a leading (A) ~ 
iV 2 behaviour, as is physically reasonable in the limit of 
large pressures, suggests that it is reasonable to assume 
that the ring should swell smoothly at all pressures, with 
the average area and surface tension described by the gen- 
eral scaling forms (A) = N 2 h{pN)^ = f 2 {pN) [7||S]. 
This argument also accounts for the possibility of a phase 
transition in the self-intersecting case, since the exponent 
of the argument in the scaling functions above diverges 
at v = 1/2. Numerical simulations of a specific general- 
isation of the LSF model support the conjecture that no 
phase transitions occur upon varying p in this model and, 
by extension, in the general case [7] . 

Consider, however, Figs, [^a) and (b). Both repre- 
sent the variation of an appropriately normalised area as 
a function of a scaled pressure variable, in two different 
models of pressurised, self-avoiding polymer rings. As is 
apparent, the average area varies discontinuously in one 
model (Fig. HJa)), whereas the other (Fig. dfb)) appears 
to show a continuous transition. 

Such behaviour is precluded by the scaling arguments 
discussed above, but it is easy to reason that the discrep- 
ancy originates in the extension of the blob picture to 
regimes where it is inapplicable. When p ~ 1, the ener- 
getics of the linkages connecting the monomers is probed 
[6]. Equivalently, the blob size reduces to the size of the 
monomer in this limit. Thus v = 1 since the ring is no 
longer fractal at or below this scale, and the exponent 
(2 - 2v)/(2v - 1) in the expansion (A) ~ N 2 (pN)^i 
vanishes. Sub-leading corrections from perimeter contri- 
butions dependent on the detailed potential between par- 
ticles constituting the ring, must then be accounted for. 
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Fig. 2: Order parameter 



for (a) [Model A], (b) [Model B] and (c) [Model C] and the scaled asphericity NA 2 for (d) [Model 



A], (e) [Model B] and (f) [Model C], illustrating their variation with the scaled pressure p. Model A, corresponding to (a) and 
(d), shows a sharp discontinuous jump in the order parameter, mirrored to a lesser extent in the asphericity. The inset to (a) 
shows the proportionality p c — k/2 for Model A. For Model B, corresponding to (b) and (e), both quantities vary continuously 
across the transition. The inset in (b) shows collapse of the data for different Rq, when p — p c and (A) are scaled as in Eq. [S] For 
Model C, as shown in (c) and (f), the transition can be sharp or non-existent depending on the ratio 8 = a/Ro- This behaviour 
can be seen both in the order parameter as well as the asphericity. 



Models. — We define and study three models of pres- 
surised self-avoiding ring polymers. The first, model A, 
consists of point particles connected by quadratic springs 
with a maximum extension of Ro. The springs cannot 
intersect, thus enforcing self-avoidance. The potential be- 
tween the particles takes the form 

V(r) = hr 2 , r<R a , 

= co, r>R . (2) 

The second, model B, again assumes point particles con- 
nected by springs with a maximum extension R$. How- 
ever, the spring potential in this case is chosen to be of the 
Finitely Extensible Nonlinear Elastic (FENE) [13] type 
and is given by 

V(r) = -Rll n (l-^j, r<R 0l 

= oo, r > Ro (3) 

This form of the potential allows us to investigate the ef- 
fect of a smooth divergence to infinity for r approaching 
Ro, as opposed to the jump imposed in model A. 

Our model C consists of monomers of finite diameter 
a connected by springs with a maximum length Rq. Self- 
avoidance implies that that no bead overlaps with another 
bead and that no bonds intersect each other. The spring 



potential is a quadratic potential with a cutoff at the max- 
imum bond length Ro, as in model A. The relevant dimen- 
sionless parameter in this model is the ratio of the bead 
size to the maximum bond- length, < 6 = < 1. In the 
limit 5 — 0, we recover model A, whereas for S = 1/1.8, we 
obtain the analog of the LSF model, albeit with springs 
instead of tethers. The spring constant is chosen to be 
k = 2; results for other values of k are obtained by rescal- 
ing the pressure axis as p — > 2p/k (see below). 

Simulation Details and Measured Quantities. 

We have performed Monte-Carlo simulations of models A, 
B and C , combining both local (single bead) and non- 
local (multiple bead) moves P3]. The basic Monte Carlo 
move consists of choosing a bead at random and giving it a 
random displacement. Self-avoidance is ensured by check- 
ing that the random displacement does not lead to bond 
intersections or bead overlaps, if applicable [15]. The sys- 
tem is evolved typically for ~ 4 x 10 6 Monte Carlo steps, 
with the relevant averages computed after discarding the 
first 10 6 steps for equilibration. We have also performed 
multi-particle moves involving the reflection of a group of 
n connected particles across the line joining the particles 
at the end. A variety of initial configurations are used to 
initiate the simulations. Defining bead positions on an un- 
derlying fine grid speeds up the simulations considerably. 

We define the order parameter as the ratio of the average 
area to the maximum allowed area, $ = A/A max with 
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Fig. 3: Susceptibility plots for (a) model A, (b) model B and (c) model C. The first-order character of the area susceptibility of 
model A is shown in (a), with a delta-function singularity at the transition. In contrast, model B displays a power-law divergence 
of the susceptibility for different values of the bond length (Ro), shown in (b). The inset to (b) shows the susceptibility of model 
B for different values of the system size, which all collapse onto a single curve. The susceptibility for model C is shown for (c) 
for an intermediate value of the bond length (Ro = 3). The plot shows two clear peaks, the second peak scaling with pN, and 
the first with pN 3 ^ 2 as shown in the inset to (c). 



A max = \NEjq cot(4r), as appropriate to an 7V-gon of 
fixed side length Rq. In the thermodynamic limit, defined 
as the double limit N — > 00 and Rq — > 00, this order 
parameter should be zero for p < p c and nonzero above it. 

The nature of the transition can be characterised 
through an appropriately defined susceptibility. We de- 
fine the area susceptibility [TTj 



X 



1 d(A) 



(4) 



For a continuous transition, both the order parameter and 
susceptibility collapse around the critical point can be de- 
scribed by appropriately defined scaling functions 



<T> = R^hx[(p-p c )R^], 
X = R^h 2 [(p-p c )R^]. 



(5) 
(6) 



The instantaneous shape of the ring is characterised by 
the gyration tensor, defined as T 0)j g = jf^2f = i(Xi. a — 
XcG, a ){Xi,p-XcG,p), with X CG , a denoting the a th com- 
ponent of the centre of mass and X^ a denoting the a th 
component of the i th particle. If the eigenvalues of this 
tensor are Ai and A2, the radius of gyration R 2 G = Ai + A2 
and the asphcricity A2 is 



An 
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Ai + A 2 



(7) 



This asphericity is for a perfectly spherical shape and 1 
for a rod-like shape. 

Results. — We first discuss results for model A. The 
order parameter shows a first order jump for the system 
sizes shown, in which the number of monomers N ranging 
from 16 to 64, as well as for two different values of the 
spring constant k. At smaller k, the transition p shifts to 
smaller values; as shown in the inset, the critical scaled 
pressure equals k/2. There are no significant finite size ef- 
fects in the data, with the data for N = 16 overlapping the 



data for N = 64; this is also the case for much larger sys- 
tems (not shown). Fig. EJd) shows the asphericity across 
the transition; while Ai vanishes asymptotically as N is 
increased both above and below the transition indicating 
that the asymptotic shape is a circle in both cases, there is 
a significant difference in the A^-dependence of A 2 on both 
sides of the transition. The susceptibility (Fig. 02[a)) dis- 
plays a classic first-order delta-function peak at the tran- 
sition (note the logarithmic scale on the y-axis) with the 
peak height increasing with N. 

The dependence of the critical pressure on k can be 
computed in the following way, relying on the assump- 
tion that fluctuations can be ignored for p > p c : The 
pressurisation energy is Fp = —4tt 2 p^-. Adding the 
spring energy Fs = -f-^T^) 2 = 27r 2 fc-^- to this gives 
Fp + Fs = [— 2p + k] and thus a first-order transition 
dX p c = k/2, as seen in the data. The self- avoidance term 
will, of course, prevent collapse for p < p c ; a Flory theory 
discussed below indicates that (A) ~ A^ 2 in this phase. 

In contrast, order parameter results for model B 
(Fig. [21(b)), indicate a smooth increase from zero at the 
critical pressure value. The order parameter curves show 
no dependence on N but depends on Rq above the tran- 
sition. The asphericity (Fig.^d)) varies strongly with N 
below the transition, while the A^-dependence essentially 
collapses above the transition. The susceptibility plots 
for model B are shown in Fig. Efb). Data for different 
values of Rq can be collapsed using Eq. [6j where the ex- 
ponents (1 « 1.33 and (2 ~ 1-33 as shown in the inset 
to Fig. [2]Jb). The order parameter variation can also be 
collapsed using the form given by Eq. [5] with exponents 
given by Q\ s» 0.67 and 82 ~ 1.33. (The errors on these 
exponents are large, however, and conservative estimates 
of error bars are around ±0.10 on each of the calculated 
exponents.) For p < p c , we have ^ ~ (p c — p)~ Bl / 92 , as 
illustrated in Fig. SJa) . Above the transition, the area or- 
der parameter can be calculated analytically (as shown in 
the next section) yielding, for p > p c , ~ (1 — l/p)R 2 , 
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Fig. 4: (a)The divergence of area with pressure near the critical 
point for a range of system sizes and maximum bond lengths 
for model B; (b) The large pressure behaviour of the average 
area for model B. The data is shown for three system sizes N 
= 128, 256, 512 with R = 30.0 . 



as illustrated in Fig. [4jb) . 

These results for self-avoiding rings with point parti- 
cles then imply the following: If there is a finite maxi- 
mum extension Rq, then there is a transition such that 
lim^Q-^co jv^oo — ls zero below the transition and non- 
zero above it. This transition separates a weakly expanded 
phase, in which the area scales as iV 2 but is otherwise 
insensitive to the value of Rq, from a strongly expanded 
phase, in which the area scales as (NRq) 2 . Depending 
on the form of the potential, the transition between these 
phases can either be discontinuous or continuous. The 
transition can also be characterised as a shape transition: 
below p, the asphericity Ai is Rq- independent for finite N 
while it decreases to zero with increasing Rq above it. 

To address the question of transitions in models of 
the LSF type, we now present our results for model C, 
where the monomers are now beads of a finite diameter a 
connected by springs with a maximum allowed extension 
Rq. Order parameter plots for this model are shown in 
Fig. [He). For the smallest Rq, the order parameter varies 
smoothly through the transition. As Rq is increased, how- 
ever, the variation becomes sharper, with what appears 
to be discontinuous behaviour of the order parameter at 
large Rq. This behaviour is also apparent in the aspheric- 
ity, shown for two values of Rq in Fig. [H[f). For Rq = 2 
(S = 1/2), there is no transition, while for i?o = 30 
(5 = 1/30), a sharp jump in the asphericity across the 



transition is seen. The susceptibility plots for this model 
(Fig. [3Jc)) show two distinct peaks, the first about p = 
scaling as pN 3 / 2 , and the second about p = 1, scaling as 
pN. 

Thus either, (a) there is a critical value of 5 below which 
there is no transition as a function of p and above which 
there is a discontinuous transition or, (b) there is no tran- 
sition except in the limit of 5 — > (model B), where a tran- 
sition demonstrably exists. We conjecture that the second 
scenario is applicable, for the following reason: As p — > p c , 
blob sizes becomes comparable to the microscopic cutoff. 
However, for finite beads, the monomer dimension bounds 
the blob size from below since fluctuations below this size 
cannot be resolved, giving the perimeter a finite width 
even in the limit of very large pressures. The smoothness 
assumption implicit in the blob-based scaling arguments 
should thus continue to hold, with both area and the sur- 
face tension varying non-singularly as p is increased. For 
models with non- vanishing 5, we therefore conjecture that 
there can only be a crossover, though possibly a sharp one, 
from the weakly to the strongly expanded phase. This is 
consistent with the observed absence of the transition at 
large pressures in LSF-likc models [7]. 

Formally, a vesicle perimeter of mean length L ~ 27ri? 
undergoes transverse thermal wandering of a magnitude 

1/2 ... 

u rms ~ (2nkBT/p) ' [6_. Note that u rrns is independent 
of N , validating a central assumption of the Flory theory 
presented later. Given 7 = pR, the ratio of the intrinsic 
width £ (blob size) to u rms is ^ ~ (pN)^ {^k^rY' 2 - 
Inserting v = 3/4, we obtain ~ (p) 3 / 2 (p/27rfc B T) 1 / 2 . 
With p -)• 0, N — ► 00, with p — Np = const, ^ma. — >• 00; 
this conclusion holds true even if we assume v — )• 1. Thus, 
any scale which bounds the blob size £ from below, such 
as a bead dimension, should provide a smooth variation 
of (A) at all pressures. 

Scaling Arguments. — A Landau formulation justi- 
fying the existence of a transition can be framed in terms 
of the competition between a surface free energy arising 
out of an effective tension term -f^^fr) an d the pressuri- 
sation term —Pjir- The expansion of the surface tension 

term yields terms of the form^(ao + 12772 + 04774 ■ • ■), 
where the lowest order contribution (ao > 0) is present 
even in the tethered case, from entropic elasticity. Com- 
bining terms of the same order yields a Landau free en- 
ergy of the form Fl ~ oiqjt- + etijb . . ., with the ao term 
changing sign at a critical value of the pressure. Such a 
Landau theory predicts an (area) order parameter expo- 
nent of 1 for the continuous case (model B), while com- 
paring free energies on both sides of the transitions yields 
first order behaviour for model A, as seen. The nature of 
the p < pc phase requires consideration of the effects of 
self- avoidance. 

A simple Flory-like theory, accounting for contributions 
from self-avoidance, pressurisation, entropy and bond 
stretching, provides a consistent explanation of the be- 
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haviour at all pressures. The stretching free energy is a 
function of the maximum bond length Rq and the form 
of the potential between neighbouring monomers. For 
model B, (and, in fact, fairly generically) it is of the form 

F stretching ~ ~fiTjy3 J,- Apart from a self-avoidance term 
FsA-, we may take over terms from the study of the pres- 
surised self-intersecting ring, to get a Flory free energy 
F ~ F S a + isr(l - P) + -RTjvJ ■ To estimate F S a, we 
assume that since the model is always in an expanded 
phase with (A) ~ A^ 2 , F$a is well represented by small 
fixed- width fluctuations about the N-gon shape. The con- 
tribution to the free energy from self-avoidance is then 
approximated as F$a = N 2 /R (as opposed to N 2 /R 2 for 
a two dimensional self avoiding ring) . Minimising this free 
energy yields the following predictions, 



i(l-p)2/3 ' V VC 

N 2 RT, p = p c (8) 

N 2 R 2 (p-l), p>p c 

Note that (A) always scales as N 2 and is independent of Rq 
below the transition. Above the transition, (A) ~ (p— 1), 
increasing linearly just above the transition, consistent 
with Fig. 2(b) and Landau-theory predictions. Thus, such 
a Flory- type scaling theory predicts the existence of a tran- 
sition at a critical value of the scaled pressure. The expo- 
nents for the area scaling predicted using this Flory theory 
are 9\ ~ 4/5 and 62 — 6/5, at the outside limit of the er- 
ror bars of the values we obtain numerically for Model B 
(0.67 and 1.33). However, errors on these exponents in 
the simulations are substantial and the predictions of our 
Flory-type theory may possibly be correct. 

The behaviour of the area above the transition can 
be computed analytically. In this regime, the shape of 
the polygon ring is simply a regular polygon and the 
free energy is simply the sum of the pressurisation en- 
ergy and the spring energy. For a regular polygon of 
side a, the area is given by A — jNa 2 cot(j^), while 
A m ax = \NR 2 cot(j^), where Rq is the maximum allowed 
bond length. For Model B, the free energy can be written 
as, F = —pA — NR 2 In (l — A A ) . The average area is 

obtained from 8F/dA = 0, giving p = 4tan (V jV ) _ Defining 
p = ipcot(-^-) , we obtain p — Wr— . Note that this 

definition for p is in fact the correct one, reducing to the 
usual p = for large system sizes. Thus 

-fi-max P 

This relation predicts the behaviour of the area very accu- 
rately for large pressures and is consistent with numerical 
data extending almost up to the critical pressure for large 
Rq , where it coincides with the prediction from Flory the- 
ory: (A) ~ N 2 R 2 (p — 1) to lowest order. This is shown in 
Fig. ffib). 



Conclusions. — We have studied the thermodynam- 
ics of several models which generalise the original model of 
LSF. We demonstrate the existence of a phase transition in 
an appropriately scaled, pressure-like quantity. This tran- 
sition is best thought of as a shape transition, in which the 
polymer conformation and area scaling are independent of 
a parameter Rq which sets the maximum bond size be- 
low the critical pressure, while depending strongly upon it 
above. The polymer ring is un-faceted below the transition 
and faceted in the form of an N-gon just above, swelling 
smoothly thereafter while retaining its faceted shape as 
the pressure is increased further above p c . Such transi- 
tions can be discontinuous or continuous in the limit that 
the ratio of the bead size to the characteristic scale of bead 
separation vanishes. In the general limit, where this ra- 
tio is non-zero, we conjecture that no transition exists, a 
conclusion which should apply to models of the LSF type. 

The possibility of a genuine phase transition in an ap- 
propriately scaled pressure variable in models for self- 
avoiding pressurised rings has not been raised before, 
although the presence of such a transition in the self- 
intersecting case has been settled decisively [TTJ[T2]. Sim- 
ple Flory-type theories appear to give useful insights into 
the transition. A similar study for the three-dimensional 
case may yield important insights and would appear to be 
very worthwhile. 

* * * 

We thank D. Dhar for useful discussions. 
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